An Estimate of the Vibrational Frequencies of Spherical 

Virus Particles 

L.H. Forcf] 

Department of Physics and Astronomy 
Tufts University, Medford, MA 02155 

Abstract 

The possible normal modes of vibration of a nearly spherical virus particle are discussed. 
Two simple models for the particle are treated, a liquid drop model and an elastic sphere 
model. Some estimates for the lowest vibrational frequency are given for each model. It is 
concluded that this frequency is likely to be of the order of a few GHz for particles with a 
radius of the order of 50nm. 

PACS numbers: 87.50.Kk, 43.80.Jz, 87.15.La 



* Email: ford@cosmos .phy.tufts.edu 



1 



I. INTRODUCTION 

Virus particles (virions) come in a variety of sizes and shapes. However, approx- 
imately spherical shapes with diameters in the range between 50nm and lOOnm are 
especially common. Many nearly spherical viruses are revealed by X-ray crystal- 
lography to have icosahedral symmetry. A typical virus particle contains genetic 
material, RNA or DNA, surrounded by a protein coat (capsid). Such an object 
should have reasonably distinct vibrational frequencies, the study of which may be 
of interest. Excitation of these vibrations could have applications in either the diag- 
nosis or treatment of viral diseases. To this author's knowledge, the sole discussion 
of these vibrational modes in the literature is that of Babincova et al [1]. These 
authors discuss the conjecture that ultrasound in the GHz range could be resonantly 
absorbed by HIV virus particles, leading to their destruction. Cooper et al Q have 
recently reported the detection of viruses by acoustic oscillations. However, the pro- 
cess of "rupture event scanning" , which these authors report, involves the separation 
of a virus particle from antibodies by ultrasound. This is distinct from the excita- 
tion of the vibrational modes of the virus particle itself, and occurs at much lower 
frequencies. 

There have also been some experimental studies of ultrasonic absorption by empty 
viral capsids 0, H| . These experiments reveal an enhanced absorption in the MHz 
range as proteins reassemble into a capsid, but do not find a resonant peak in this 
frequency range. Witz and Brown |j| have emphasized that these and other results 
show that viral capsids are flexible and change size or shape in response to vibrations 
or to changes in temperature or pH. 

The purpose of the present paper is to provide some estimates of the lowest vibra- 
tional frequencies of a spherical virus particle. The simplest estimate is to take this 
frequency to be of the order of a characteristic speed of sound divided by the size of 
the virus particle. This is the estimate used in Ref. For the purpose of giving a 
more accurate estimate, we will examine two models, which treat the particle (1) as 
a liquid drop and (2) as a uniform elastic sphere. Similar models have been used by 
Bulatov et al [6] to estimate the vibrational frequencies of nanoclusters. 



II. A LIQUID DROP MODEL 

Consider a sphere of radius a filled with a nonviscous liquid with surface tension 7 
and mass density p. The lowest vibrational mode of this sphere will be a quadrupole 



mode with frequency 0] 





which can be written as 
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Liquid 


/ / 1 vv 


01 ' OVJ 
r 1 r vv 


v/(10 s Hz) 


Benzene 


0.397 


0.88 


2.3 


Diethylene glycol 


0.62 


1.12 


2.5 


Trehalose 


0.95 


1.63 


2.6 


Lysine hydrochloride 


0.90 


1.38 


2.7 


Arginine hydrochloride 


0.95 


1.44 


2.8 



TABLE I: The mass density, surface tension, and the lowest vibrational frequency predicted 
by Eq. © for drops of various liquids with a radius of a = 50nm. The data for Benzene 



and Diethylene glycol 



are for droplets in air at room temperature. The data for the 



three proteins |9j are for aqueous solutions at approximately 50° C. 



where pw = 10 3 kg/m 3 and jw = 0.073Nm are the mass density and surface tension 
for water, respectively. The surface tension and mass density, along with the lowest 
vibrational frequency derived from Eq. (j2J) for a = 50nm, are given in Table 1 for 
several liquids. 

Recall that Eq. (j2J) assumes a nonviscous liquid. In fact, the viscosity of many 
of the liquids in Table 1 cannot be neglected for such a small drop. When viscosity 
is sufficiently large, the drop will not oscillate, but rather undergo overdamped mo- 
tion [7J. The main lesson from the liquid drop model is that a nonviscous liquid drop 
of a = 50nm with a typical surface tension and mass density would have a lowest 
vibrational frequency of the order of a few times 10 8 Hz. 



III. AN ELASTIC SPHERE MODEL 

A better model for a virus particle is to treat it as a uniform elastic sphere. The 
three independent parameters which characterize such a sphere can be taken to be 
the radius a, the speed of pressure waves, cp, and the speed of shear waves, cs- The 
oscillations of an elastic sphere are treated in detail by Pao and Mow Here we 
quote their results in the notation of Ref. jgj. The eigenfrequencies of the normal 
modes are given by the vanishing of the determinant 

S rp (n,x,y) S rs (n,y) _ 

— U , (6) 

S tp (n,x) S ts (n,y) 

where 

S rp (n, x, y) = (n 2 -n - ^y 2 ) 2 j n (x) + 2xj n+1 {x) , (4) 
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FIG. 1: The root x of Eq. © is plotted as a function of b = cg/cp for n = (compression 
mode) and n = 2 (quadrupole deformation). 

S TS (n,y) = n(n + l)[(n - - xj n+ i(y)] , (5) 

S tp (n, x) = (n- l)j n (x) - xj n+ i(x) , (6) 

S ts (n,y) = -(n 2 - n - ^x 2 ) 2 j n (x) - xj n+1 {x) . (7) 
Here x = ua/cp, y = ua/cs, and the j n are spherical Bessel functions. Let 

6=- (8) 

be the ratio of the speed of the shear wave to that of the pressure wave, so that y = bx. 
Given b and n, we can solve Eq. © for x and hence for the frequency of the normal 
mode, v = xcp/(2ira). The smallest roots for x for the n = and n = 2 modes are 
plotted in Fig. d The corresponding frequencies of oscillation can be expressed as 

, = 4.8xl0 9 ^f^ ( ° P , )x. (9) 

The frequencies obtained from this equation for the n = mode are given in Table 
II for various materials. 

We can see that for a wide range of materials, the lowest mode of vibration is a 
purely radial mode with a frequency of the order of a few times 10 9 Hz for a sphere 
of radius a = 50nm. 
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Material 


Cp 


cs 


v/(W 9 Hz) 

I \ / 


Nylon [111 


2620 


1070 


11.4 


Polystyrene [111 


2350 


1120 


11.3 


Polyethylene [11] 


1850 


540 


5.6 


Neoprene rubber [llj 


1600 




5 x 


Polynucleotides [12] 


1700 - 1900 




(5.3 - 6) x 


Amino acids [12] 


1900 - 2400 




(6 - 7.7) x 


Globular proteins [12] 


1700 - 1800 




(5.3 - 5.7) x 



TABLE II: The speeds of sound (in m/s) and the lowest vibrational frequency for the n = 
mode predicted by Eq. @ for spheres of various materials with a radius of a = 50nm. In 
cases where no data for cs are available, the frequency is given as a multiple of x, which is 
likely to be less than one for these materials. 

IV. CONCLUSIONS 

In the previous sections, we have examined two models for a spherical virus par- 
ticle, a liquid drop model and an elastic sphere model. It is of interest that the two 
models yield estimates for the lowest vibrational frequency which differ by only about 
one order of magnitude. Of the two models, the elastic sphere model is probably the 
better description of a virus particle. An even better model might be one in which 
the particle has a liquid core (DNA or RNA) surrounded by an elastic outer shell 
(the capsid). Such a model would probably yield vibrational frequencies intermediate 
between those predicted by the two models discussed in this paper. In any case, we 
obtain an estimate for the lowest vibrational frequency of the same order of magni- 
tude as that given in Ref. [H, in the range of a few GHz for particles with a size of 
about lOOnm. Of course, the existence of a resonance requires that damping be below 
the critical value above which overdamped motion occurs. Even if this condition is 
fulfilled, it is difficult to predict the width of the resonance. This remains a question 
for experimental investigation. The existence of well defined resonances could prove 
valuable both for basic science and for medicine. Thus this is a potentially fruitful 
area for further research. 
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